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ABSTRACT 



Spectral fitting of the spin a = cJ/GM 2 in the microquasar GRS 1915 + 105 estimate values higher than a = 0.98. However, there are 
certain doubts about this (nearly) extremal number. Confirming a high value of a > 0.9 would have significant concequences for the 
theory of high-frequency quasiperiodic oscillations (HF QPOs). Here we discuss its possible implications assuming several commonly 
used orbital models of 3:2 HF QPOs. We show that the estimate of a > 0.9 is almost inconsistent with two hot-spot (relativistic 
precession and tidal disruption) models and the warped disc resonance model. In contrast, we demonstrate that the epicyclic resonance 
and discoseismic models assuming the c- and g- modes are favoured. We extend our discussion to another two microquasars that 
display the 3:2 HF QPOs. The frequencies of these QPOs scale roughly inversely to the microquasar masses, and the differences in 
the individual spins, such as a = 0.9 compared to a = 0.7, represent a generic problem for most of the discussed geodesic 3:2 QPO 
models. To explain the observations of all the three microquasars by one unique mechanism, the models would have to accommodate 
very large non-geodesic corrections. 
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1. Introduction 

In the past few years an impressive amount of work has been 
done on estimating the black hole spin using the X-ray contin- 
uum fitting method (see, e.g., McClintock et al., 2006, 2008; 
Middleton et al., 2006; Done et al., 2007; Shafee et al., 2008; 
McClintock et al., 2010, 2011). The obtained results indicate 
that the spins of the individual sources cover almost the whole 
range spanning from a non-rotating Schwarzschild black hole to 
near-extreme Kerr black holes (McClintock, 2010; McClintock 
et al., 2010, 2011). The latter type of object has been detected 
in the microquasar GRS 1915+105 for which it was estimated 
by McClintock et al. (2006) that the dimensionless spin a = 
cJ/GM 2 is higher than a = 0.98. However, some doubts remain 
about the (nearly) extremal value reported by McClintock et al. 
In particular, using a similar fitting procedure Middleton et al. 
(2006) derived a value of a — 0.7. Recent independent spin esti- 
mates of GRS 1915+105 based on another, so called, relativistic 
iron line profile fitting imply the spin either around 0.6 or 0.98 
(Blum et al., 2009). There are several concerns about both of the 
spectral (continuum and iron-line) spin estimate methods. It is 
not the purpose of our paper to list and analyze these reasons in 
detail but it is useful to recall some findings of a direct relevance 
to observations of low-mass X-ray binaries (LMXBs). 

Within continuum method studies, the X-ray spectrum is as- 
sumed to be the sum of black body contributions from different 
disc radii. Within the adopted concept, the expected spectrum is 
determined by the radial distribution of disc temperature, which 
follows from the relativistic version of the Shakura-Sunyaev ac- 
cretion disc model and depends on the black hole spin (Shakura 
& Sunayev, 1973; Novikov & Thorne, 1973). The spin is then 
inferred by fitting the observed disc spectra to those calculated. 
In more detail, one of the key theoretical assumptions of this 



method lies in the expectation that the inner edge of the ac- 
cretion disc coincides with the innermost stable circular orbit 
(ISCO). Consequently, the spectral fitting provides an estimate 
of the ISCO position that depends on the black hole spin (see 
Shafee et al., 2006; McClintock & Remillard, 2009; Penna et 
al., 2010, for a more complex picture). 

The method deals only with a subset of observations where 
the assumption of the blackbody emission seems valid. Thus, 
only datasets that show (presumably) weakly Comptonized 
spectra acquired for the thermal dominant source state have been 
used in several works (e.g., McClintock & Remillard, 2003, 
2009; Shafee et al., 2008, but see also Steiner et al. (2009)). 
Moreover, only observations where the source luminosity is less 
than ~ 20% of the Eddington luminosity provide reliable results, 
since at high luminosities the standard ISCO concept is not valid 
(see Done & Davis, 2008; Abramowicz et al., 2010). The method 
is also sensitive to the accuracy of the source distance, mass, and 
inclination determined from independent methods (usually from 
the optical measurements). 

The iron-line method studies are based on the evidence of a 
broad, skewed iron-line profile that is also observed in the spec- 
tra of several LMXBs. The observed iron line is believed to orig- 
inate from the reflection of hard X-ray photons in the innermost 
parts of the accretion disk. These photons, emitted by an external 
process (often assumed to be inverse-Compton scattering in a hot 
corona), have sufficient energy to remove the K-shell electrons 
from the iron atoms in the accreted gas and produce the Fe-K 
emission. The broad profile of the Fe-K line is then explained 
by relativistic Doppler effects, light bending, and the gravita- 
tional redshift characteristic of the vicinity of the binary central 
compact object (see Miller, 2007, for a review). As for spectral 
continuum fitting, it is assumed that the disc extends down to 
ISCO. Consequently, properties of the broadened emission line 
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are associated with the ISCO radii when considering the spin as 
a free parameter. The advantage of this method is that it does not 
require accurate information about the source distance and mass, 
and provides estimate of the source inclination itself (e.g., Davis 
et al., 2006; McClintock & Remillard, 2009). Its weak points are 
connected namely to important requirements about the strength 
of the detected signal and to difficulties in subtracting the line 
from the whole spectra (when adopting detailed modeling of the 
continuum). Moreover, there are doubts that the line broadening 
is caused by general relativistic effects, since these strongly de- 
pend not only on the radius of the emission, but also on details 
of the local emission that are not very known (e.g., Beckwith & 
Done, 2004). 

Both of the aforementioned means of estimating the black 
hole spin rely on being able to relate the source states to the disc 
geometry. In the case of GRS 1915+105, it is unclear how this 
can be achieved (e.g., McClintock & Remillard, 2003; Done et 
al., 2004; van Oers et al., 2010; Li et al., 2010). In view of all 
the uncertainties, one should keep in mind that some caution is 
needed in relation to the high value of the black hole spin in 
GRS 1915+105. 1 

2. Orbital models of quasiperiodic oscillations 

The X-ray power density spectra (PDS) of several LMXBs ob- 
tained during the past three decades contain peaked features 
called quasi-periodic oscillations (QPOs). These QPOs are ob- 
served in neutron star (NS) as well as black hole (BH) systems 
and occur with periods of the order of 1 - 10~ 3 sec. An outstand- 
ing example of these oscillations arises in the high frequency 
part of the PDS being thus called "high frequency" (HF) QPOs. 
When more HF QPOs are observed in black hole sources, their 
frequencies seem to form commensurable pairs with a preferred 
ratio of 3:2 (Abramowicz & Kluzniak, 2001; McClintock & 
Remillard, 2003). The black hole HF QPOs are typically weaker 
and less coherent than most of the NS HF QPOs (e.g., Barret 
et al., 2005a,b; Mendez, 2006). The neutron star HF QPOs ap- 
pear as two correlated modes at distinct (but, in contrast to 
BH HF QPOs, variable) frequencies. The range of frequencies 
spanned by the two NS modes is rather large (ranging from 50Hz 
to 1500Hz, being typically a few hundreds of Hz for each indi- 
vidual source). The two NS modes are in a way analogous to the 
BH case because they usually exchange their dominance when 
passing the 3:2 frequency ratio (Torok et al., 2008a,b,c; Torok, 
2009; Boutelier et al., 2010). Properties of the HF QPOs are re- 
viewed in, e.g., McClintock & Remillard (2003) or van der Klis 
(2006). 

There is strong evidence that the HF QPOs originate very 
close (less than 100 gravitational radii r g = GMc~ 2 ) to the accret- 
ing compact objects, but no commonly accepted QPO theory has 
so far been developed (e.g., van der Klis, 2006). Moreover, while 
there are several concrete, well-described similarities and differ- 
ences between the black hole and neutron star HF QPO phe- 
nomenology (including those briefly mentioned above), it has 
not yet been resolved whether the generic mechanism could be 
the same for both classes of the sources. There are several pro- 
posed models that cannot be applied to both the classes, e.g., 

1 McClintock et al. (2011) states that a > 0.98 but the authors also 
note that some caution is needed namely because of the incomplete 
study of the source distance (see their paper for details). In this pa- 
per, we show that confirming it would have a significant impact on 
the theory of the high-frequency quasiperiodic oscillations (HF QPOs) 
and present the discussion of some concrete implications for a few fre- 
quently quoted QPO models. 



due to the requirement of a solid NS surface (Lamb et al., 1985; 
Alpar & Shaham, 1985) or a high black hole spin (Stuchlfk et al., 
2007a,b; Slany & Stuchlfk, 2008). Nevertheless, a larger variety 
of models has been designed under the assumption of a common 
mechanism (see, e.g., van der Klis, 2006, for a basic overview). 
Most of the hypotheses assume a relation between the QPO fre- 
quencies and the frequencies related to motion of accreted matter 
orbiting in the vicinity of a compact object (hereafter "orbital" 
models). These models typically deal with either "hot-spot" or 
"disc-oscillation" QPO interpretation. We next consider an arbi- 
trary choice of several commonly quoted models. For the sake 
of comprehensibility, we first give a short summary of the exam- 
ined models and recall some of their main features and related 
references. 

2. 1 . Kinema tic models 

The kinematics of the orbital motion allows the consideration 
of the variability that arises from the motion of "hot-spots" or- 
biting inside the accretion disc. We consider two models of this 
kind. The "relativistic precession" model (hereafter RP model) 
was proposed in a series of papers by Stella & Vietri (1998a, 
1999, 2002); Morsink & Stella (1999). The model illustrates 
that the kHz QPOs represent modes of the relativistic epicyclic 
motion of blobs at various radii r in the inner parts of the ac- 
cretion disc. It is often recalled and known for roughly match- 
ing the correlation between the HF QPOs observed in the NS 
sources (e.g., Belloni et al., 2007a; Torok et al., 2010; Lin et al., 
2010). Within the model, the twin-peak QPO frequency correla- 
tion arises because of the periastron precession of the relativis- 
tic orbits. Owing to Lense-Thirring relativistic precession, the 
model also predicts another frequency correlation extending to 
higher timescales. The kHz QPO frequencies are indeed corre- 
lated to the low-frequency QPO features observed around ~1- 
50Hz (e.g., Stella & Vietri, 1998b, this correlation is however 
not being the focus of our present work). On the other hand, 
there are apparent difficulties with the detailed modeling of the 
relation between the hot-spot motion and the observed modula- 
tion (e.g., Lamb & Markovic, 2000). 

Cadez et al. (2008), Kostic et al. (2009), and Germana et al. 
(2009) introduced a similar concept in which the QPOs are gen- 
erated by a "tidal disruption" (TD) of large accreting inhomo- 
geneities. Their hydrodynamic simulations signify that blobs or- 
biting the central compact object can be stretched by tidal forces 
forming "ring-section" features within the model expected to be 
responsible for the observed modulation. The model has been 
proposed for both supermassive and stellar mass black holes . We 
note that, at least in some cases, the power density spectra simu- 
lated using the model closely reproduce those that are observed. 
However, the concept of "rocks" approaching the radii of the or- 
der of ISCO-radius in LMXBs is questionable. The upper limit 
to the radii, where the tidal forces begin to disrupt the object of 
density p orbiting a primary source of mass M, is roughly de- 
scribed by the Roche-limit rpD ~ (M/p) 1 ^ 3 . The related charac- 
teristic orbital frequency vtd ~ (GM/r^) 1 ' 2 can be expressed as 
vtd = (Gp) 1 ^ 2 . For the rocks it is clearly given by vtd ~ 10 3 Hz, 
which is far lower than the observed (~ 100Hz) frequencies. 2 

Identification of the lower and upper kHz QPO frequencies 
with the frequencies of the orbital motion is recalled in the top 
part of Table 1 for both hot spot models. To match the observed 
3:2 ratio, the QPOs in the RP model must be generated at the 
radii where v K /v r = 3/1 (r = 6.75GMc~ 2 for a = 0), i.e., very 

2 We thank the anonymous referee for emphasizing this issue. 
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close to ISCO (r = 6GMc 2 for a — 0) where the radial epicyclic 
frequency vanishes. For the TD model, this location is shifted 
outwards close to the radius where the radial epicyclic frequency 
exhibits a maximum value (v r /v r = 2/1; r = 8GMc 2 for a = 0). 

2.2. Resonant models 

Within the (present) RP model, there is no generic explanation 
of the observed 3:2 frequency ratio R = 3/2 and there is a 
clear need to explore the issue of the preference of certain or- 
bits. The concept of the TD model is less problematic because 
for any a and r it implies that R e (1,2) and the predicted ef- 
fects occur around the location of the maximal allowed v r , where 
v r /v r ~ 2/1. Nevertheless, the frequency ratio is also not re- 
liably constrained by this model. In contrast, the commensu- 
rability of the frequencies is crucial for the models when as- 
suming the warped disc (WD) oscillations suggested by Kato 
(2001, 2004a,b, 2005, 2008). Kato (2008) reviewed several pos- 
sible resonances in deformed disc (see also works of Ferreira & 
Ogilvie, 2008, 2009). The main weak point of the concrete WD 
model assumed here (Kato, 2004a,b, see Table 1) is that it con- 
siders a somewhat exotic disc geometry that causes a doubling 
of the observed lower QPO frequency. 

The frequency commensurability is also crucial for the non- 
linear resonance models discussed by Abramowicz, Kluzniak 
and collaborators (Kluzniak & Abramowicz, 2000; Abramowicz 
& Kluzniak, 2001; Abramowicz et al., 2003a,b; Bursa et al., 
2004; Rebusco, 2004; Horak, 2008; Stuchlik et al., 2008a,b; 
Horak et al., 2009, and others; see also Aliev & Galtsov, 1981 
and Aliev, 2007). A particular resonance of this kind that is of- 
ten discussed is the 3:2 internal epicyclic or Keplerian resonance 
(Ep, Kp, see Torok et al., 2005, for details). The explicit formu- 
lae of the frequency relations corresponding to these three reso- 
nant disc oscillation models are listed in the bottom part of Table 
1. We note that while for the WD resonance model the QPOs 
are located at the same radii as for the kinematic TD model, 
i.e., around the location of the maximum of v r , the two 3:2 reso- 
nances occur above this maximum, at a larger distance from the 
central object. We also consider here another two QPO (reso- 
nance) models that we denote as RP1 model (Bursa, 2005) and 
RP2 model (Torok et al., 2010). Both of them assume different 
combinations of non-axisymmetric disc-oscillation modes. They 
are of particular interest because they involve oscillation modes 
whose frequencies for slow rotation almost coincide with the 
frequencies predicted by the RP model. These two resonances 
should occur much closer to ISCO than the resonance expected 
in the Ep model. 

The physical interpretation of the RP1 model is unclear, 
since the oscillation modes assumed within the model are un- 
likely to be able to enter the resonance (see, Horak, 2008, for 
details). The resonant coupling between the pairs of the oscilla- 
tion modes assumed in the Kp, EP, and RP2 models is in prin- 
ciple allowed, but detailed physical mechanisms providing this 
coupling and mode excitations have not yet been fully devel- 
oped (see, e.g., Sramkova et al., 2007; Horak, 2008; Rebusco, 
2008; Kluzniak, 2008, and references therein). Substantial effort 
clearly yet needs to be invested in extending both the related an- 
alytic and numerical work and the link between them. 

2.3. Models assuming fundamental discoseismic modes 

We have so far only considered models in which it is assumed 
that both of the observed 3:2 frequencies are produced by the 



same mechanism and excited at a certain (common) preferred 
radius. A qualitatively different consideration relating each of 
the two frequencies to a different radial region arises when gen- 
eral combinations of the fundamental discoseismic modes are 
assumed. The HF QPOs have been proposed to correspond to the 
following three distinct modes: the so-called g-modes (inertial- 
gravity waves that occur at the radius where the radial epicyclic 
frequency reaches its maximum value), c-modes (corrugation 
vertically incompressible waves near the inner edge of the disk), 
and p-modes (inertial-pressure oscillations that occur near the 
edge of the disc). This concept was elaborated in particular by 
Kato & Fukue (1980), Okazaki et al. (1987), Nowak & Wagoner 
(1992), Wagoner (1999), Silbergleit et al. (2001), Wagoner et al. 
(2001), Ortega-Rodriguez et al. (2002), Silbergleit & Wagoner 
(2008), and Wagoner (2008). Disc oscillation modes that have 
been proposed to explain the QPOs have been observed in hy- 
drodynamical (HD) simulations of the accretion processes (e.g., 
Zanotti et al., 2005; Reynolds & Coleman, 2009) but there is 
growing evidence from magnetohydrodynamic (MHD) simula- 
tions that these modes are typicaly dumped by various instabil- 
ities produced in the presence of the magnetic field (see, e.g., 
Tsang & Lai, 2009; Fu & Lai, 2009, 201 1). At the same time, the 
MHD simulations do not convincingly reproduce the 3:2 QPOs 
that undoubtedly appear in the X-ray fluxes of the LMXBs. It is 
therefore also unclear whether these simulations accomodate all 
the crucial ingredients or not (see, e.g., Machida & Matsumoto, 
2007, 2008, who found in their MHD simulations some low fre- 
quency QPOs only in the presence of the radiative cooling). 

2.4. Confronting the QPO models with the rapid spin of black 
hole in the microquasar GRS 1915+105 

In general, the existing QPO models represent rather unfinished 
concepts that have specific advantages and difficulties. The os- 
cillations predicted by these models are usually not seen in 
the present MHD simulations. Moreover, none of these models 
yet really match the full QPO phenomenology observed in the 
LMXBs including the dependence of the QPO visibility on the 
source spectral states, QPO amplitudes, and coherence times. In 
this situation, we expect that the implications of the eventually 
confirmed high spin of GRS 1915+105 could help us improve 
the individual QPO models or even discriminate between them. 
In the next section, we calculate the spin implied by the above- 
mentioned QPO models. 

3. The spin implied by individual models 

There is a straightforward connection between the hot-spot mod- 
els and the characteristic orbital frequencies (i.e., the three fre- 
quencies of the perturbed circular orbital motion: the azimuthal 
"Keplerian" frequency, the radial epicyclic frequency, and the 
vertical epicyclic frequency). The frequencies of numerous disc- 
oscillation modes are also constrained by the three characteristic 
orbital frequencies. Assuming a Kerr geometry, the orbital fre- 
quencies for a given radius depend only on mass and spin of the 
black hole. 3 It is therefore possible to infer the black hole spin or 
mass from the observed 3:2 frequencies and (at least some) con- 
crete orbital models. This procedure has been previously per- 
formed, e.g., for a wide selection of the resonance QPO mod- 
els (Abramowicz & Kluzniak, 2001; Torok et al., 2005; Torok, 

3 We consider here only the Kerr spacetime as a standard description 
for rotating BHs, although alternatives have been discussed in a similar 
context (see Kotrlova et al., 2008; Stuchlik & Kotiiova, 2009). 
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Fig. 1. Ratio of the frequencies of the fundamental discoseismic modes (based on Wagoner et al., 2001). Left: Ratio of v c to v g . The spread of 
two functions corresponds to the uncertainty in the speed of sound. The horizontal dotted line denotes the 3:2 frequency ratio. Values in the boxes 
evaluate the spin ranges required by this ratio. Right: The same consideration, but including the p-modes. 

Table 1. Frequency relations corresponding to individual QPO models and the spins implied by the 3:2 QPOs in GRS 1915+105 and the mass- 
range 10- 18Mq. The middle column indicates the ratio of the epicyclic frequencies determining the radii corresponding to the observed 3:2 ratio. 
The indicated ranges of spin also represent total spin ranges for the whole group of the three microquasars (see Figure 2). 
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2005). We recall that the complete set of formulae required to 
evaluate the Keplerian and epicyclic orbital frequencies in Kerr 
geometry was first derived by Aliev & Galtsov (1981). These 
frequencies were extensively discussed in several later studies. 
A detailed analysis can be found in Torok & Stuchlfk (2005), 
where the definition formulae are given in the usual form and 
the frequencies are studied across the full range of positive a. 
We note that here we do not consider a > 1 since we focus our 
attention fully on models originally designed for BHs and not on 
their possible extensions to naked singularities. 

The 3:2 QPO frequencies in GRS 1915+105 are well known 
to be given by (McClintock & Remillard, 2003) 

vu = 168 + 3Hzandv L = 113 + 5Hz. (1) 

Assuming Eq. (1) and the aforementioned formulae for the 
orbital frequencies, we calculate the implied mass-spin func- 
tions for the models associating the 3:2 QPOs with a com- 
mon radii by means of the definition relations given in Table 1 . 
Following Abramowicz & Kluzniak (2001) and Torok et al. 
(2005) and taking into account the estimated range of the mass 
of GRS 1915+105 (e.g., McClintock & Remillard, 2003) 

10 M <M< 18 M , (2) 

we infer the expected ranges of the spin. We found that the RP, 
TD, WD and RP2 models imply that the spin is rather too low, 
i.e. a < 0.6. Somewhat more satisfactory are the predictions ob- 
tained from the Ep (0.65 < a < 1), Kep (0.7 < a < 1), and RP1 



(a < 0.8) models. For each model, the detailed results are pre- 
sented in the last column of Table 1 . As recalled above, for the 
discoseismic modes the individual observed QPOs correspond 
to different modes located at their own radii. The frequencies 
of these modes depend on the black hole spin and the speed of 
sound in the accreted gas, and scale roughly as 1 /M, whereas 
their dependence on the other parameters of the accreting sys- 
tem is supposed to be very weak (Wagoner et al., 2001). The 
frequencies of the c- and g- modes (v c and v g ) depend only mod- 
erately on the speed of sound. Thus, their ratio is mainly a func- 
tion of the black hole spin, which is depicted in the left panel of 
Figure 1 . The right panel of Figure 1 shows the relations of the 
ratios of the frequencies v c or v g to the p- mode frequency v p , 
which are strongly dependent on the speed of sound. Inspecting 
both panels of Figure 1, we can see that for the combination of 
modes relating v v = v p and v L = v a the frequency ratio R is never 
higher than unity. For the other five possible combinations of 
modes, there is in each case some range of a that corresponds 
to the 3:2 frequency ratio. This range is located above a = 0.9 
only for the combination relating v v = v c and v L = v g . In more 
detail, this combination requires that 0.9 < a < 0.94. 4 The mass 
needed to match v v = 168 + 3Hz is then 13M Q < M < 19M Q 



4 For higher values of a, the model predictions seem to deviate from 
the 3:2 ratio since the c-mode frequencies, when a increases, grow faster 
than the frequencies of the g-mode. However, the published works have 
not fully investigated the behaviour for very high spin values a > 0.95 
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with good overlap 13M 
given in Eq. (2). 



< M < 18M with the expected mass 



4. Discussion and conclusions 

It follows from our previous discussion that the hypothesis of 
the internal (epicyclic) resonance 5 and the discoseismic model 
which relates the upper and lower 3:2 QPO to the c- and g- 
modes are favoured in the case of GRS 1915+105 provided 
that a > 0.9. The TD, WD, RP, and RP2 models are then dis- 
favoured. This statement was inferred assuming that v K , v r , and vg 
are the exact geodesic frequencies. A similar analysis including 
the influence of non-geodesic effects would require very detailed 
study. Here we only roughly estimate the possible relevance of 
the non-geodesic effects. Figure 2 shows the v u xM{a) curves im- 
plied by the individual geodesic models. These curves (as well 
as the curves drawn in the consequent Figure 3) have non-trivial 
ambiguous extensions for the range of a > 1 representing the 
naked-singularity region. We do not consider them here, since 
we focus our attention on black holes (we expect to discuss the 
subject of a > 1 in a different work). The observationally de- 
termined interval of v v x M for GRS 1915+105 is indicated by 
the light yellow rectangle. We define the relative non-geodesic 
correction to be 



— (^observed ^predicted) /^predicted? 



(3) 



which is needed to match the observations of GRS 1915+105 
with a given model for a certain spin. From Figure 2, we can 
easily estimate by eye that the non-geodesic correction Av w re- 
quired for the frequency predicted by the RP model for a = 0.95 
is about -50%. From Figure 3, we can find that for a e (0.9, 1), 
the quantity Av w changes from -40% to -60%. The same is 
roughly true for the TD and WD models, while for the RP2 
model the required correction is even higher. Thus, our result 
is justified, except when considering of very large non-geodesic 
corrections. 
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Fig. 2. Curves v L ,xM(a) implied by the individual geodesic models. The 
light yellow rectangle indicates the observationaly determined interval 
of v v x M for GRS 1915+105. The red dotted vertical line denotes a = 
0.95. The red vertical arrow indicates the correction needed to match the 
upper limit to v v x M with the RP model for this spin. The colour boxes 
are drawn for the mass and spectral spin estimates given by different 
authors for GRS 1915+105, GRO J1655-40 and XTE 1550-564 (see 
the main text for references). The dotted blue line indicates the lower 
observational limit to Mx v v that is roughly common to GRO J1655-405 
and XTE 1550-564. 



4.1. Relation to other sources 

Torok et al. (2005) pointed out that since the 3:2 QPO frequen- 
cies in microquasars scale roughly as v v =2.8(M/M Q ) _1 kHz 
(McClintock & Remillard, 2003), their spins implied by the 
epicyclic resonance model should not vary much among them. In 
Figure 2, we include several examples of spin estimates obtained 
by different groups and methods for GRS 1915 + 105 and the 
other two microquasars, GRO J 1655-40 and XTE 1550-564. The 
references related to these estimates are as follows (the assumed 
ranges of v v being taken from McClintock & Remillard, 2003). 
For the mass, Greene et al. (2001), Greiner et al. (2001), Orosz 
et al. (2002), Beer & Podsiadlowski (2002), and McClintock & 
Remillard (2003) provide commonly accepted mass estimates. 
These are in Figure 2 denoted by letter A. Beer & Podsiadlowski 
(2002) present an alternative prediction (in Figure 2 denoted by 
letter B) that moves the lower boundary of estimated mass of 
GRO J1655-40 from 6.0 to 5.1 M . For the spin, there were the 



so far because of numerical difficulties (Wagoner 2010, from a private 
communication) . 

5 Torok et al. (2005) assumed this and several other resonances for 
the spin estimates of the Galactic microquasars. From the other consid- 
ered resonances, the "5:1" resonance can also match the high spin. This 
resonance involves relatively high resonant coefficients and occurs very 
close to ISCO. 



following studies: [1] McClintock et al. (2006); [2] Middleton et 
al. (2006); [3] McClintock et al. (2008); [4] Miller et al. (2009). 

The above quoted spin estimates assume either the spectral 
continuum or the iron line method. It is apparent from Figure 2 
that both different authors and different methods suggest (some- 
what) different values of spin (and for one microquasar even dif- 
ferent values of mass). The epicyclic resonance model favoured 
(along with the discoseismic model) in the GRS 1915+105 
seems to match at least some of these estimates, but for the pa- 
rameters of the XTE J1550-564 (a -0.7, M~5 - 1M Q ) assumed 
in the figure it fails. If different spins (a > 0.9 in GRS 1915+105 
and a ~ 0.7 in GRO J1655-40 and XTE J1550-564) were con- 
firmed, the difficulty of matching the all observed 3:2 frequen- 
cies would clearly be rather generic for most of the orbital QPO 
models. 

The observationally determined ranges of v v x M for GRO 
J1655-40 and XTE J1550-564 nearly coincide, the upper limits 
being roughly equal to the upper limit for GRS 1915 + 105 (see 
Figure 2). Taking advantage of this setup, we plot a simple rough 
scheme of the non-geodesic corrections required for a given 
model, spin, and source in Figure 3. The red curves in the figure 
indicate the minimal corrections required for GRS 1915 + 105. 
Parts of these curves with a negative sign also roughly indicate 
the corrections required for the other two microquasars. Positive 
corrections for these microquasars are then included as the blue 
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curves. If the values a > 0.9 for GRS 1915+105 and a ~ 0.7 
for GRO J 1655-40 or XTE J 1550-564 were confirmed, the re- 
quested corrections would appear to be rather high for most of 
the denoted models assuming a unified 3:2 QPO mechanism. 
Only the RP1 model can survive with corrections of |Av| up to 
~ 20% (but as recalled in Section 2, the present physical inter- 
pretation of this model is questionable). 

Because of the observational 1/M scaling, the above diffi- 
culty also arises for the discoseismic models (which are not con- 
sidered in Figures 2 and 3). For these, we present in Table 2 
the mass ranges implied by all combinations of the fundamental 
modes. These appear to overlap well with those observationally 
determined only for the model relating the upper and lower 3:2 
QPO to the c- and g- mode provided that a e (0.90 - 0.94). For 
the other combinations related to different spins, the mass ranges 
differ from those in the observation. Clearly, there is the need for 
a substantial correction also for a unified 3:2 QPO model assum- 
ing fundamental discoseismic modes provided that microquasars 
had different spins of a > 0.9 compared to a ~ 0.7. 
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